Modeling of Disk-Star Interaction: Different Regimes of Accretion and Variability 

Marina M. Romanova, Akshay K. Kulkarni 
Dept. of Astronomy, Cornell University, Ithaca, NY 14853 
akshaySastro . Cornell . edu , romanovaOastro . Cornell . edu 

Min Long 

Dept. of Physics, Center for Theoretical Astrophysics, University of Illinois at Urbana-Champaign, 

UrbanalL, 61801 

longSastro . Cornell . edu 

Richard V.E. Lovelace 

Dept. of Astronomy and Dept. of Applied and Eng. Physics, Cornell University, Ithaca, NY 14853 

lovelaceOastro . Cornell . edu 

ABSTRACT 

The appearance and time variabiHty of accreting miUisecond X-ray pulsars (hereafter 
AMXPs, e.g. Wijnands & van der Klis 1998) depends strongly on the accretion rate, 
the effective viscosity and the effective magnetic diffusivity of the disk-magnetosphere 
boundary. The accretion rate is the main parameter which determines the location of 
the magnetospheric radius of the star for a given stellar magnetic field. We introduce 
a classification of accreting neutron stars as a function of the accretion rate and show 
the corresponding stages obtained from our global 3D magnetohydrodynamic (MHD) 
simulations and from our axisymmetric MHD simulations. We discuss the expected 
variability features in each stage of accretion, both periodic and quasi-periodic (QPOs) . 
We conclude that the periodicity may be suppressed at both very high and very low 
accretion rates. In addition the periodicity may disappear when ordered funnel flow 
accretion is replaced by disordered accretion through the interchange instability. 

Subject headings: accretion, accretion discs; instabilities; MHD; stars: oscillations; stars: 
magnetic fields 



1. Different Regimes of Accretion 



A. BOUNDARY LAYER (BL). If the accretion rate is sufficiently large that Vm < R*, then the star's 
magnetic field is completely buried by the accreting matter (e.g. Gumming, Zweibel, Bildsten 2001; 
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Fig. 1. — Main regimes of accretion to a rotating neutron star which occur for different accretion 
rates. Cases A to F show different regimes when the accretion rate gradually decreases. A is the 
Boundary Layer regime where the magnetic field is completely buried or screened, and matter 
accretes directly to the surface of the star; B is the Magnetic Boundary Layer regime where the 
magnetosphere is small but may influence the djmamics of the matter flow; Regimes C and D are 
where the magnetosphere is large enough to channel the matter flow; however, accretion may also 
occur through instabilities (case C); Regime E is the Propeller stage. Regime F is the radio-pulsar 
stage. 
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Lovelace et al. 2005), and matter accretes to the star directly through a boundary layer, which 
is expected to be approximately axisymmetric because there is no channeling by the magnetic 
field. In this case, high-frequency QPOs may be associated with inhomogeneities at the boundary 
between inner radius of the disk and the star. Axisymmetric hydrodynamic simulations of accretion 
through the BL have been recently performed by Fisker & Balsara (2005). Additional factors such as 
radiative pressure in the inner regions of the disk around NS may also determine the characteristic 
radius in the disk and corresponding high-frequency QPOs (e.g.. Miller, Lamb, & Psaltis 1998; Miller 
&Lamb 2001). 

B. MAGNETIC BOUNDARY LAYER (MBL). At lower accretion rates, the star's magnetic field is only 
partially buried, so that it influences matter flow around the star. We will draw an approximate 
line between MBL case and case of larger magnetospheres at = 2i2*, so that the MBL cases 
correspond to < Vm < 2i?*. Even in this range there are a number of possible cases with 
different observational properties (Romanova & Kulkarni 2008). 

MBL with No Magnetosphere. If the magnetospheric radius is very close to the stellar radius, then 
it cannot channel matter to the poles. However, some magnetic flux threads the inner regions 
of the disk and thus may influence the dynamics of the disk (e.g., Warner & Woudt 2002). The 
authors suggested that the inner region of the disk may form an equatorial belt, which may be 
responsible for the QPOs in dwarf novae, a sub-class of cataclysmic variables. Such MBLs were 
proposed theoretically and should be investigated in greater detail in the future. First exploratory 
simulations have been done by Romanova & Kulkarni (2008) . 

MBL with Small Magnetosphere. We investigated accretion to a star with a small magnetosphere 
in global 3D MHD simulations (Romanova & Kulkarni 2008). We observed from simulations that 
matter may accrete in stable or unstable regime. 

First we investigated accretion in stable regime to a star with misalignment angle 9 = 15° 
in two cases of small and tiny magnetosphere. We observed that small magnetosphere formed 




Fig. 2. — 3D MHD simulations of accretion through the magnetic boundary layer of a star (in 
the stable regime) with a misalignment angle Q = 15°. Left two panels show slices of density 
distribution and sample magnetic field lines in the case of a small magnetosphere. The two right- 
hand panels show similar slices but for a tiny magnetosphere. 
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and matter accreted through funnel streams (see Figure 2, left panels) . In the case of about twice 
as smaller dipole field of the star a tiny magnetosphere formed and also disrupted the disk (see 
Figure 2, right panels) . Equatorial slice shows that the magnetic field lines are wrapped by the disk 
and may influence its structure. In particular, non-axisymmetric bending-type waves have been 
observed in the inner parts of the disk. 

Another run at smaller misalignment angle = 5° led to unstable accretion but of unusual 
kind: instead of sporadically forming unstable tongues as in the unstable regime (Kulkarni & Ro- 
manova 2008a,b; Romanova, Kulkarni & Lovelace 2008), we observed the formation of two sym- 
metric equatorial tongues which connect the star with the disk and which rotate around the star 
with angular velocity of the inner disk (see Figure 3). The tongues form two equatorial spots which 
rotate faster than the star: each spot rotates about 4 times faster than the star. The frequency 
approximately corresponds to the frequency of the inner disk and will change when position of 
the inner disk change. Thus a star in this accretion regime will reveal high frequency QPOs with 
frequency drift determined by the accretion rate to the star. Such frequency drift and correlation 
with matter flux is observed in AMXPs (e.g., van der Klis 2000). 

C and D. UNSTABLE AND STABLE REGIMES OF ACCRETION. For even lower accretion rates, the 
disk matter is stopped by the magnetosphere at a larger distance, a few radii of the star. Recent 3D 
MHD simulations have shown that in this case the matter may accrete either through symmetric 
funnel streams which are lifted above the magnetosphere (stable regime), or through an inter- 
change instability, where temporary stochastic tongues of matter form at the disk-magnetosphere 
boundary and penetrate deeply into magnetosphere. These tongues form multiple hot spots on the 
surface of the star. Figure 4 shows both regimes of accretion and the corresponding light-curves 
associated with the hot spots on the surface of the star. 

Stable Regime of Accretion (Funnel Streams). In the stable regime, matter accretes in two ordered 
funnel streams (Figure 4, left panel, see 3D MHD simulations by Romanova et al. 2003, 2004a; 
Kulkarni & Romanova 2005). These two funnels form two ordered spots at the surface of the 
star, which vary only slightly around their preferred positions, so that the stellar period is clearly 




Fig. 3. — 3D MHD simulations of accretion through an unusual magnetic boundary layer which 
formed in the case of a small misalignment angle 6 = 5°. From left to right: 3D image of the 
accretion disk; an equatorial slice of the density distribution, and the Fourier power spectrum. 
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Type 


M imMQ/yr) 


Type of Variability 




Boundary Layer 


M > 7.3 X 10-8 


QPO 




Magnetic BL 


6.5 X 10^^ < M < 7.3 X 10^8 


QPO (PJ 


2i?* < < 3.1i?* 


Funnel, Instabilities 


1.4 X 10"^ < M < 6.5 X 10"^ 


P*, QPO 


3.1i?* < r„ < 12R, 


"Propeller" 


1.3 X 10^11 < M < 1.4 X 10-9 


QPO (PJ 


> 12R^ 


Pulsar 


M < 1.3 X 10-^^ 


radio P^, 



Table 1 : Ranges of magnetospheric radii, accretion rates and expected variability features in differ- 
ent regimes of accretion for a neutron star with P* = 2.5ms. 




Fig. 4. — Stable (left) and unstable (right) regimes of accretion. The bottom panels show light- 
curves from the hot spots (from Romanova, Kulkarni & Lovelace 2008) . 
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observed in the Fourier power spectrum. This is the regime where X-ray pulsations in AMXPs 
are expected. If the magnetic field is more complex than the dipole, then the hot spots have a 
more complex structure. However, spots will have a definite preferred locations so that a definite 
periodicity is expected (Long et al. 2008). Our simulations shown that a star may spin up or 
down due to interaction of the magnetosphere with the disk, so that the star's period is expected to 
wander around the rotational equilibrium state, which corresponds to the condition ~ l-2rcor, 
which has been derived from simulations by Long et al. (2005) (see also Ghosh & Lamb 1979). 
Such variations of the period have been recently observed, e.g., by Burderi et al. (2007) and Papitto 
et al. (2008). In addition to periodic variations associated with the hot spots, QPOs associated with 
oscillations of the inner disk regions are expected (Lovelace & Romanova 2007; Lai & Zhang 2008). 

Simulations have shown that at small misalignment angles 9 < 5° a hot spot may rotate faster 
or slower than the star depending on the accretion rate (Romanova et al. 2003; Romanova et al. 
2006). This may lead to the rapid phase variations and other interesting phenomena observed in 
AMXPs (Lamb et al. 2008a,b). 

The Unstable Regime of Accretion has been observed in multiple 3D MHD simulations in cases where 
the misalignment angle of the dipole is small, 6 = 5° (Kulkami & Romanova 2008a,b; Romanova et 
al. 2008). In the unstable regime, matter accretes to the star through several "tongues" which form 
at the inner edge of the disk and accrete to the surface of the star on the dynamical time-scale, 
forming temporary hot spots on its surface. The light-curve from the hot spots looks stochastic 
(see Figure 4, right panel), and in many cases no period is detected in the wavelet or Fourier 
spectra. Thus, in the strongly unstable regime, a star with a full-sized magnetosphere may be in 
the stochastic regime of accretion. In less strongly unstable cases, both tongues and funnels are 
present, and periodicity can be detected. 

Boundary Between Stable and Unstable Regimes: We observed from numerous simulations that (1) 
the stable regime usually dominates at high misalignment angles, Q > 30°; (2) enhancement of 
the accretion rate usually leads to the onset of the unstable regime (Kulkarni & Romanova, these 
Proceedings). The last condition can be interpreted as follows: for a NS with a given period of 
rotation P*, the growth of the accretion rate leads the gravitational force to dominate over the 
centrifugal (or, << ?'cor)> and the instability appears more easily. At smaller accretion rates, the 
magnetospheric radius moves out, closer to rcor, and the centrifugal force becomes strong enough 
to oppose instability. Detailed comparisons of stable/unstable regimes with instability criteria were 
also performed (Kulkarni & Romanova 2008a; Romanova et al. 2008). From other side we observed 
that at larger accretion rate the density of matter at the inner edge of the disk is larger which is 
favorable for the onset of instability. 

F. "PROPELLER" REGIME. At even smaller accretion rates, the magnetospheric radius rm becomes 
larger than the corotation radius rcor, and the star enters the "propeller" regime. Axisymmetric 
simulations performed in our group have shown that propeller regime may be either strong or 
weak. 
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Fig. 5. — Disk-magnetosphere interaction for different magnetic Prandtl numbers. The left-hand 
panel is for Pr = 1, when the disk matter flows in and penetrates through the magnetosphere 
diffusively at the same rate. The right-hand panel is for Pr >> 1, where the radial accretion speed 
is significantly larger than the diffusion speed of the field through the matter. 




Fig. 6. — The left panel shows outflow of matter (background, selected field lines (red) and stream- 
lines (blue) in the strong propeller regime. The fluxes at the right panel show that the process of 
accretion and outflows is episodic (from Romanova et al. 2005). 
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In the Strong Propeller Regime, disk matter acquires angular momentum from the rotating magne- 
tosphere fast enough that most of it is ejected in conical outflows. At the same time a significant 
amount of angular momentum and energy flows along the open stellar field lines, giving an axially- 
directed Poynting-flux dominated jet (Romanova et al. 2005; Ustyugova et al. 2006). Multiple 
simulations have shown that the strong propeller operates when Pr » 1. In this case the field 
lines are bunched into an X-type point (see Figure 5, right panel) which is favorable for outflows, as 
discussed by Shu and collaborators (e.g., Cai et al. 2008). In this regime, the cycles of matter accu- 
mulation at the inner disk, its diffusion through magnetosphere and outbursts to jets are observed 
with some accretion to the star (see also Goodson et al. 1997). Smaller ratios r„Jr^.„,. and higher 
diffusivities (still at Pr 1) are also favorable for the strong propeller regime. Oscillations are 
often quasi-periodic, with frequency vqpq = 0.01 — 0.2v^, where the lowest frequency is observed 
at the lowest diffusivity, and vice-versa. 

In the Weak Propeller Regime, when Pr 1, the magnetosphere interacts with the disk, and the star 
loses its angular momentum without outflows. The disk may oscillate due to such interaction, and 
most of the incoming disk matter may accrete to the star's surface (Romanova et al. 2004b). It is 
often the case that the rotating magnetosphere (modified by the azimuthal wrapping of the field) 
pushes the disk to larger distances, and matter may accrete from that point. In the limit of a very 
weak propeller where ~ rear, stellar pulsations and different QPO frequencies are expected as in 
the case of funnel flow accretion. If magnetosphere is large enough then accretion may be stopped 
by the centrifugal barrier of the fast rotating magnetosphere in both "propeller" regimes. 

G. PULSAR REGIME. For even smaller accretion rates, the radius of the magnetosphere becomes 
larger than the light cylinder radius, = c/il^ . In this regime, accretion to the surface of the NS 
is suppressed, and the star becomes a radio-pulsar. This stage is least investigated numerically and 
it is not clear, for example, at what values of Vm/rLc the radio-pulsar emission dominates. 



2. Discussion 

Analysis of different stages of evolution of accreting magnetized neutron stars have shown that 
they have different appearances for different accretion rates. In addition, their appearance may be 
different depending on other factors, such as the magnetic diffusivity and the Prandtl number 
(ratio of viscosity to diffusivity). We chose as an example a typical NS with period = 2.5 ms 
{v^ = 400Hz), and derived accretion rates from equation (2) corresponding to the boundaries 
between different regimes. The results are shown in the Table. As we see from simulations of 
different regimes, periodic oscillations (with the star's period) are observed only in the regime of 
magnetospheric accretion which corresponds to a rather narrow range of accretion rates, 1.4 x 
10~^ Mq/pv < M < 6.5 X 10~^Mo/yr. This fact may explain the appearance of pulsations in 
intermittent AMXPs in a narrow range of M as discussed by Altamirano et al. (2008) in the case of 
SAX J1748. 9-2021. If the pulsations are also present in the magnetic boundary layer case and in 
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the weak propeller regime, then these boundaries for M will be somewhat wider. 
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